Azimuthally fluctuating magnetic field and its impacts on observables in heavy-ion collisions 
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The heavy-ion collisions can produce extremely strong transient magnetic and electric fields. We study the 
azimuthal fluctuation of these fields and their correlations with the also fluctuating matter geometry (character- 
ized by the participant plane harmonics) using event-by-event simulations. A sizable suppression of the angular 
correlations between the magnetic field and the 2nd and 4th harmonic participant planes is found in very cen- 
tral and very peripheral collisions, while the magnitudes of these correlations peak around impact parameter 
b ~ 8 — lOfm for RHIC collisions. This can lead to notable impacts on a number of observables related to 
various magnetic field induced effects, and our finding suggests that the optimal event class for measuring them 
should be that corresponding to b ~ 8 — 10 fm. 

PACS numbers: 12.38.Mh, 25.75.-q, 25.75.Ag, 24.10.Nz 



I. INTRODUCTION 

Ultra-relativistic heavy-ion collisions create not only a do- 
main of extremely high energy density where a new state 
of matter — the deconfined quark-gluon plasma (QGP) may 
form, but also extremely strong (electro)magnetic fields due 
to the relativistic motion of the colliding heavy ions carry- 
ing large positive electric charge [1]. Previous computations 
showed that the magnetic fields generated in Au + Au colli- 
sion at RHIC (s/s = 200GeV) can reach about eB ~ m\ ~ 
10 18 G [2-8], which is 10 13 times larger than the strongest 
man-made steady magnetic field in the laboratory. The mag- 
netic field generated at LHC energy can be an order of mag- 
nitude larger than that at RHIC [3, 8], according to a simple 
scaling law recently found in [8] for the event-averaged mag- 
netic field: (eB y ) oc Zby/s for b < 2Ra , where Z is the 
charge number of the ions, b is the impact parameter, Ra is 
the radius of nucleus, and the j/-axis is perpendicular to the 
reaction plane. Thus, heavy-ion collisions provide a unique 
terrestrial environment with ultra-strong magnetic fields. 

There have been very strong interests and intensive efforts 
recently in studying various possible physical effects induced 
by the presence of strong magnetic field. Receiving partic- 
ular enthusiasm is the set of ideas to look for experimental 
manifestation of QCD effects stemming from topology and 
anomaly and aided by the external magnetic field. These in- 
clude, e.g., the so-called Chiral Magnetic Effect (CME) [2, 9- 
11] in which a nonzero axial charge density in the matter 
(presumably from topological transitions via sphalerons) to- 
gether with the magnetic field B will induce a dipole charge 
separation along the B direction. A lot of works have been 
done in the past few years to experimentally measure this ef- 
fect by analyzing the charged particle correlations [12-14] 
and to understand the interpretation and background effects 
related to these measurements [15-20], as recently reviewed 
in Ref. [21]. There is also the so-called Chiral Separation Ef- 
fect (CSE) [22], from its interplay with the CME there arises a 
gapless collective excitation called the Chiral Magnetic Wave 
(CMW) [23]. An observable effect of CMW was proposed in 
[24]: the CMW with the presence of nonzero vector charge 



density and the magnetic field transports the charges in QGP 
towards an electric quadrupole distribution with more posi- 
tive charges near the poles of the produced fireball (pointing 
outside of the reaction plane) while more negative charges 
near the equator (in the reaction plane), and this leads to a 
measurable splitting of negative and positive pions' elliptic 
flow. Recently STAR collaboration reported the first measure- 
ment of the charged pion flow spitting versus the charge asym- 
metry which appears in agreement with the predictions from 
CMW [24, 25]. Yet one more interesting effect is possible 
soft photon production through the QCD conformal anomaly 
in the external magnetic field as suggested in [26]. In this 
mechanism the photons are emitted perpendicular to the B di- 
rection and there is an appreciable azimuthal anisotropy of the 
emitted photons which might partially account for the unusu- 
ally large v-i of direct photons reported by PHENIX collabo- 
ration [27]. In addition to the above anomaly phenomena in 
magnetic field, there are also discussions of other novel ef- 
fects in strong magnetic fields [28-31], e.g., the spontaneous 
electromagnetic superconductivity of QCD vacuum, the pos- 
sible enhancement of elliptic flow of charged particles, pos- 
sible anisotropic electromagnetic radiation from the QGP, the 
energy loss due to the synchrotron radiation of quarks, and the 
emergence of anisotropic viscosities. 

In most of the phenomenological studies of these effects, 
the magnetic field B in heavy-ion collisions represents the 
largest source of uncertainty in their quantitative calculations, 
and a precise knowledge of the magnetic field is urgently 
needed. This is particularly so in the context of the realiza- 
tion and intensive investigations in the last two or three years 
that there are very strong fluctuations in the initial conditions 
of heavy-ion collisions. Such fluctuations have been shown 
to lead to strong observable effects both in the bulk collective 
expansions (as "harmonic flows") [32] and in the anisotropy 
of penetrating hard probe (as "harmonic tomography") [33]. 
Since the magnetic field directly relies on the initial distribu- 
tions of protons in both nuclei, a realistic computation has to 
take into account such strong initial fluctuations. A first step 
has been made in [7, 8] to compute the magnitude of B with 
event-by-event fluctuations, which was indeed found to be re- 
markably modified from the naive "optical" estimates. 
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In this paper, we focus on the strong fluctuations in the az- 
imuthal orientation of the magnetic field B, and particularly 
investigate its angular correlation with the underlying mat- 
ter geometry (specified by participant planes) bearing con- 
current fluctuations on an event-by-event basis. This corre- 
lation is a very important link between the experimental mea- 
surements and any of the magnetic field induced effects in 
heavy-ion collisions. Most (if not all) previous studies rely 
on the assumption that the magnetic field is pointing in the 
out-of-plane direction which would be true without fluctua- 
tions. In the real world, however, both the B and the "planes" 
(event planes or participant planes) bear strong fluctuations 
and as we will show their orientations are never fully "locked" 
and only strongly correlated in some circumstances. We will 
study these correlations in great details and examine the con- 
sequences for a number of observables related to some afore- 
mentioned effects. 

The rest of this paper is organized as follows. In Sec. II 
we show how the initial fluctuations modify the magnetic and 
electric fields generation in heavy-ion collisions. We study 
the azimuthal fluctuation of the magnetic field in Sec. III. We 
discuss the physical implications of these results in Sec. IV. 
For completeness we also show results for the electric field 
in Sec. V. In sec. VI a new class of charge dependent mea- 
surements partly motivated by our results will be discussed. 
Finally we summarize in Sec. VII. The natural unit H = c = 
ks = 1 will be used throughout this article. 

II. EVENT-BY-EVENT CALCULATION OF THE 
ELECTROMAGNETIC FIELD 

The aim of this section is to give an event-by-event calcula- 
tion of the electromagnetic fields in heavy-ion collisions. We 
focus on the fields at the initial time t = 0, that is, the moment 
when the two colliding nuclei overlap completely. Our start- 
ing point is the Lorentz boosted Coulomb formulas which are 
equivalent to the Lienard-Wiechert potentials for constantly 
moving charges: 

eB( *' r) = £E^(Rn) K _ ( ^"^ )2]3/2 V»xR ni 

eE ^ - ^E^^) K -(R,::t„ )2] 3/ 2 R - ^ 

where R„ = r r n (t) is the relative position of the field 
point r to the nth proton at time t, r n (t), and v n is the ve- 
locity of the rith proton. The summations run over all protons 
in the projectile and target nuclei. Equations (2.1)contain sin- 
gularities at R n = if we treat protons as point charges. In 
practical calculation, to avoid such singularities we treat pro- 
tons as uniformly charged spheres with radius R p . The charge 
number factor Z n (R a ) in Eqs. (2.1) is introduced to encode 
this aspect: when the field point locates outside the nth proton 
(in the rest frame of the proton) Z n = 1, otherwise Z n < 1 
depends on R n . The in-medium charge radius R p of proton is 
unknown (the most recent measurement of the rms charge ra- 
dius of proton gives R p = 0.84184(67) fm in vacuum [34]), 
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FIG. 1 : (Color online) The event-averag ed (eB) 2 and (eE) 2 (in unit 
of m*) at t = and four different points on the transverse plane as 
functions of the impact parameter b. 

we choose R p = 0.7 fm in our numerical simulations. We 
checked that varying R p from 0.6 fm to 0.9 fm will shift the 
numerical results within 15% but no qualitative conclusion is 
altered. The nucleons in one nucleus move at constant ve- 
locity along the beam direction (we choose it as z-direction), 
with the nucleons in another nucleus moving oppositely at the 
same speed. The energy for each nucleon is set to be y/s/2 in 
the center-of-mass frame, therefore the value of the velocity 
of each nucleon is given by vf t = 1 — {2m^ / ^/s) 2 , where 
7njv is the mass of the nucleon. We set the x-axis along the 
impact parameter vector so that the reaction plane is the x-z 
plane. Finally, the positions of nucleons in the rest frame of a 
nucleus are sampled according to the Woods-Saxon distribu- 
tion. 

In this computation we have obtained in each event all 
the components of both the B and E fields at four different 
points, r = (0,0,0),r = (3fm,0,0),r = (0,3fm,0), and 
r = (3 fm, 3 fm, 0), on the transverse plane for a wide range 
of impact parameter b in Au + Au collision at RHIC energy 
\/s = 200 GeV. We have found that all the z components are 
negligibly small, while the transverse components are strong, 
with their centrality trends in agreement with previous results 
in [7, 8]. The absolute magnitudes of these fields are some- 
what smaller than those reported in [7, 8], due to the impor- 
tant difference that we treat each proton as a charged sphere 
with radius R p rather than a point charge, which is both more 
physical and mathematically less of singularity problem. As 
also noted in [7, 8] and confirmed in our computation, al- 
though the x-component of the magnetic field as well as the 
x- and y-components of the electric field vanish after aver- 
aging over many events, their magnitudes in each event can 
be huge and comparable to the 7/-component of the magnetic 
field due to the fluctuations. This already implies that both 
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FIG. 2: (Color online) The event-by-event histograms of — ^2 at impact parameters b = 0, 5, 10, 12 fm for Au + Au collision at RHIC 
energy. Here 'I'b is the azimuthal direction of B field (at t — and r = (0, 0, 0)) and \&2 is the second harmonic participant plane. 
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FIG. 3: The scatter plots on ^3-^2 plane at impact parameters b = 0, 5, 10, 12 fm for Au + Au collision at RHIC energy. Here \I/b is the 
azimuthal direction of B field (at t — and r = (0, 0, 0)) and $ 2 is the second harmonic participant plane. 



the magnitude and the direction of the electromagnetic field 
(albeit in the unmeasurable x-y plane) fluctuate strongly. In 
Fig. 1 we show ((eB) 2 ) and ((eE) 2 ) as functions of b at dif- 
ferent points: clearly the electromagnetic fields bear consid- 
erable inhomogeneity on the transverse plane. We also notice 
that the electric field can be very strong particularly for more 
central collisions. The ((eB) 2 ) is particularly interesting be- 
cause the signal strengths of several magnetic-field-induced 
effects are proportional to (eB) 2 , as we will discuss more in 
Sec. IV. 



III. AZIMUTHAL CORRELATIONS BETWEEN 
MAGNETIC FIELD AND MATTER GEOMETRY 

As already mentioned before, on the event-by-event basis 
the electromagnetic field fluctuates strongly both in magni- 
tude and in azimuthal direction. The direction of the B-field 
is very important as the B-field induced effects inherit this 
information and occur either along (e.g. CME, CMW) or per- 
pendicular (e.g. photon emission) to B. 

Even more important from the measurement perspective 
is the question of the B orientation with respect to a frame 
that could be identified experimentally rather than to the ideal 
world reaction plane. In reality what could be determined is 
the final-state hadrons' distribution geometry in momentum 
space, in particular the second harmonic V2 event plane (EP) 
which is relatively more tightly correlated with the initial mat- 
ter distribution geometry that can be specified primarily by the 
second harmonic £2 participant plane (PP). Therefore the re- 
ally useful information is the B orientation with respect to the 
matter geometry per event. Here we make a first detailed ex- 
amination of such a kind, studying the azimuthal correlation 



between the magnetic field and the participant planes. 

In order to do that, we determine the participant planes 
for various harmonics from the Monte Carlo Glauber simu- 
lations of the initial condition and analyze the angular cor- 
relations between the B and the participant plane orienta- 
tions from the same event. The nth harmonic participant 
plane angle \&„ and eccentricity e n are calculated from par- 
ticipant density p(r) as in the literature (e.g. [32]): eie 1 * 1 = 
-(/ epT ±P (rx)f±e i *)/(J dV^K), and e„e m *" = 
-(jrd 2 r ± p(r_ L )r"e in ^)/(/d 2 r_ L p(r_ L )r™) for n > 1. In 
this Section we focus on the 2nd harmonic participant plane 
^2 as it is the most prominent anisotropy from both geometry 
and fluctuations. Correlations of B with other harmonics will 
be discussed in Sec. IV. By determining B and ^2 in each 
event we can examine the distribution of their relative angle 
over many events: in Fig. 2 we plot the event-by-event his- 
tograms of — *2 at b = 0, 5, 10 and 12 fm for Au + Au 
collisions at y/s = 200 GeV where is the azimuthal direc- 
tion of the magnetic field at t = and r = (0,0,0). Strikingly 
at b = the histogram of $b — ^2 is almost uniform indicat- 
ing that ^b an d ^2 are basically uncorrected. For b = 5, 10 
and 12 fm, the histograms have Gaussian shapes peaking at 
7r/2 with the corresponding widths not small at all. The width 
first decreases when b increased from 5 to 10 fm and then in- 
creases again toward b = 12 fm. So although the B field 
indeed points at 7r/2 with respect to 9 2 plane on average, it 
fluctuates significantly in each event. The correlation between 
^b and ^2 is the strongest in middle-centrality while weak- 
ens much in the most central and most peripheral collisions. 

To further reveal the correlation pattern between ^b and 
'J 2 and to understand better the non-monotonous centrality 
trend of the widths of the histograms in Fig. 2, we show the 
scatter plots from all events at given b on the ^3-^2 plane 
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in Fig. 3 which visualize the 2D probability distribution den- 
sity. Again for 6 = 0, the events are almost uniformly dis- 
tributed indicating negligible correlation between and ^/ 2 . 
For b = 5, 10, and 12 fm, the event distributions evidently 
concentrate around (^^,^2) = (t/2, 0) indicating a cor- 
relation between the two. Going from b = 5 to 10 and to 
12 fm, the spread in ^b keeps shrinking while the spread in 
if? 2 clear grows with larger b. This is because for non-central 
collisions with increasing 6, the B is mostly from the specta- 
tors whose number increases and bears less fluctuations while 
the ^2 is determined by participants whose number decreases 
and fluctuates more. This explains the non-monotonic trend 
of the widths in the histograms in Fig. 2. In short, we have 
found that the event-by-event fluctuations of the initial condi- 
tion bring azimuthal fluctuations in both "Pb and *f?2, and the 
angular correlation between them is smeared out significantly 
in the very central and very peripheral collisions while stays 
strong for middle-centrality collisions. This observation cer- 
tainly influences the interpretation of observables related with 
B-induced effects, as will be discussed in the next section. 



IV. IMPACT ON VARIOUS OBSERVABLES 

We now discuss the impacts of the azimuthal fluctuations 
of the magnetic field with respect to matter geometry (the par- 
ticipant planes) on a number of pertinent observables recently 
measured in heavy-ion collisions. We recap a few points al- 
ready discussed: the B-induced effects are sensitive to the az- 
imuthal direction of the B field, which however is not exper- 
imentally known; in the past when observables are proposed 
and interpreted for measuring certain B-related effects, it is 
often assumed that the B direction is perpendicular to the re- 
action plane; as already shown in last section, this assumption 
is not true, and the azimuthal orientation between ^b of B 
and ^2 fluctuates with sizable spread in their relative angle 
(^b — ^2)- In what follows we evaluate the impacts of such 
fluctuations on the experimentally measured quantities. 

First, let us consider the pair correlation 7 = (cos(</>i + 
<f>2 — 2^2)) with 0i 2 the azimuthal angles of the particle 1 
and 2 where the average is taken over events. The measure- 
ments of 7 are motivated by the search for the Chiral Mag- 
netic Effect [2]. (The 7 is actually extracted through three 
particle correlations (cos(0i + 4>2 — 2^3)) divided by elliptic 
flow V2 with the third particle serving as a "projector" toward 
the Event Plane [12]. ) Now let us focus on the same-charge 
pairs, and suppose the CME indeed gives rise to the same- 
side azimuthal correlations along the B-field direction for the 
same-charge pairs. The two-particle density then receives the 
following contribution (with A ++ the signal strength) 

/++ = A ++ cos(0i - * B ) cos(0 2 - * B ). (4.1) 

This translates into the following form after re-defining angles 

4>i = fa - ^ 2 and <&b = *b - *2: 



We therefore see that the CME's contribution to the 74 
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If the B-direction were to be always perfectly aligned with 
the out-of -plane, i.e. ^b = tt/2, then we simply have 
7 ++ — » —(A ++ )/2. But the fluctuations in magnetic field 
as well as in matter geometry will blur the angular rela- 
tion between the two and modify the signal by the factor 
~ (cos(2^'b))- Similarly, if one measures the charge sepa- 
ration with respect to higher harmonic participant plane, e.g. 
the fourth harmonic plane ^4, the corresponding correlation 
will be (cos[2(0i + <p2 — 2^4)]) (as recently proposed with 
the hope to disentangle the collective-flow and CME con- 
tributions to 7 [35]). The azimuthal fluctuation of B field 
with respect to *f?4 will again contribute a modification fac- 
tor ~ (cos[4(5'b — ^4)]) m the above correlation. 

Recently it has been suggested that a splitting between the 
elliptic flow of the ir~ and tt + could occur in heavy-ion colli- 
sions due to the Chiral Magnetic Wave [24, 25], as briefly dis- 
cussed in Sec. I. Specifically the CMW will induce an electric 
quadruple in the net charge distribution along the B field: 

Pe {4>) ~ 2r e cos[2(0 - B )] = 2r e cos[2(0 - * B )], (4.4) 

where r e oc / dd>p e (d>) cos[2(0 — ^b)] quantifies the quadru- 
ple from CMW. This modifies the final charged hadron distri- 
bution by an amount proportional to r e : 

dN±(<f>) „ r_ . _ ,_ T , , r, 



N ± {l + 2v 2 cos(20) ± ^ cos[2(^ - * B )]} , (4.5) 



where N = N + + iV_ is the charge multiplicity. Thus one 
obtains a splitting of the charged pions V2 



-Vn = - 



cos(2tf B )>- 



(4.6) 



Again, we see the modification factor ~ (cos(2^E'b)) arising 
from relative orientation between B and ^2- If were to be 
simply 7r/2 (i.e. B always out-of-plane) then the sign of split- 
ting is always > v 2 + - However due to the fluctuation, 
deviates from 7r/2 and leads to a smaller splitting. 
Finally, we consider the example of recently suggested soft 
photon production from conformal anomaly in the presence 
of an external B field [26]. The photon emitted via this mech- 
anism carries azimuthal information of the external magnetic 
field. If one assumes a homogeneous external B field, then the 
produced photon spectrum goes like q dN/d 3 q ~(qx B) 2 
leading to the following azimuthal distribution: 

dN-y/ddi oc {1 — cos[2(0 - * B )]} = {1 - cos[2(> - # B )]} (4.7) 

which will contribute the following to the elliptic flow of these 
photons (with Vy being the signal strength) 



-(V 7 cos(2# B )>. 



(4.8) 



Once again we see the same modification factor ~ 
(cos(2^b)) due to the mismatch between B field and the out- 
of-plane direction from event to event. If ^b were to be sim- 
ply 7r/2 (i.e. B always out-of-plane) then indeed these pho- 
tons will have positive and sizable elliptic flow. However the 
fluctuation may bring ^b to depart from ir/2 considerably. 
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FIG. 4: (Color online) The correlations {cos[n(5'B — ^»)]) as functions of impact parameter for n = 1, 2, 3, 4 at four different positions on 
the transverse plane: (from left to right) r = (0, 0, 0) fm; r = (3, 0, 0) fm; r = (0, 3, 0) fm; r = (3, 3, 0) fm. 
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FIG. 5: (Color online) The (eB) -weighted correlations 
((eB) 2 cos[n(* B - *„)])/((eB) 2 ) as functions of impact 
parameter for n = 1, 2, 3, 4 at r = (0, 0, 0) fm. 

In all three examples, the fluctuating \I>b — &2 brings in a 
reduction to the intrinsic strength of the signal by 



Ri = (cos(2* B )). 



(4.9) 



It is therefore worth a close examination of this factor: see 
Fig. 4 for the computed average values of (cos[n(^B — ^n)]) 
for varied centralities from event-by-event determination of 
the B-field direction "Pb (at several different spatial points) 
and the participants harmonics, n = 1,2, 3, 4. The plots 
suggest that the correlations between ^b and the odd harmon- 
ics Vfi, ^3 are practically zero (in accord with parity invari- 
ance), while the the correlations of ^b with even harmon- 
ics ^2, ^4 are nonzero but get suppressed as compared with 
the results in optical geometry limit (without fluctuations) e.g. 
(cos[2(^b — v &2)])opt = — 1. The centrality dependence of 
(cos[2('3>b — ^2)]) is in perfect agreement with the patterns 
seen in the histograms Fig. 2 and scatter plots Fig. 3: it is sig- 
nificantly suppressed in the most central and most peripheral 
cases (indicating no correlations) while is maximized around 
b = 8 ~ 10 fm with peak values —0.6 ~ 0.7. Similar be- 
havior is observed for (cos[4(<I'b — ^4)]) too albeit with a 
weaker correlation strengths, e.g. with the peak values ~ 0.2 
or so. Although suppressed, the correlation between \I>b and 
^4 appears still sizable for moderate values of b. This may 
imply complications and caveats that must be seriously ad- 
dressed for the proposal in Ref. [35] to use the two-particle 
correlation with respect to the fourth harmonic event plane to 
disentangle the CME contribution and the flow contribution. 
The plots in Fig. 4 for four different field points demon- 



strate the spatial dependence of the correlations (cos^^e — 
$„)]) on the transverse plane. We first see minor differences 
when the field point deviates from the origin and such differ- 
ences become sizable for the outer most point r = (3, 3,0)fm. 
The comparison suggests that further away from the origin the 
correlations between ^b and $ n become even weaker. 

The last issue we want to address is the possible correlation 
between the B orientation and its own strength. It is conceiv- 
able that there could be correlation (on an event-by-event ba- 
sis) between the signal strength and the B field direction. Typ- 
ically the signal strength of these B-field induced effects will 
scale as ~ B 2 . For example, recalling that the CME current is 
proportional to B and the signal strength A ++ as mentioned 
earlier would be nearly proportional to B 2 . The quadrupole 
strength r e through CMW also scales similarly r e oc B 2 . So 
does the signal strength V 1 in the case of soft photon emis- 
sions. If the field strength and the orientation of B do cor- 
relate with each other, then one may not be able to factorize 
the signal strength and the field orientation factor cos(2^b) 
when taking the event average. We therefore further examine 
the correlations between the two by evaluating the factor 



Ri 



((eB) 2 cos(2§ B )) 
((eB) 2 ) 



(4.10) 



In Fig. 5 we show the impact parameter dependence 
of the (eB) 2 -weighted correlations ((eB) 2 cos[ti(5'b — 
*„)])/((eB) 2 ) for n = 1,2,3,4. We find little difference 
between the two factors R\ in Fig.4(the most left panel) and 
i?2 in Fig. 5 at all centralities for all n values. While only re- 
sults for field point r = (0, 0, 0) are shown in Fig. 5, we have 
checked all other three points and the observation is the same. 
Therefore, we conclude that the magnitude of the magnetic 
field has no noticeable correlation to its azimuthal direction 
with respect to the matter geometry. 



V. AZIMUTHAL CORRELATION BETWEEN ELECTRIC 
FIELD AND MATTER GEOMETRY 

In this section, we briefly discuss how the event-by-event 
fluctuation affects the correlation between the orientation of 
E-field and the participants harmonic planes, & n . The moti- 
vation is that from Fig. 1 we notice that the electric field can 
be as strong as the magnetic one, and thus they may lead to 
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FIG. 6: (Color online) The correlations (cos[n(\&E — ^n)]) as functions of impact parameter for n = 1, 2, 3, 4 at four different positions on 
the transverse plane: (from left to right) r = (0, 0, 0) fm; r = (3, 0, 0) fm; r = (0, 3, 0) fm; r = (3, 3, 0) fm. 



observable effects too. An obvious example is possible mul- 
tiple charge distributions induced by the strong electric field. 
Similarly to the magnetic field case, possible E-induced ef- 
fects would be affected by the azimuthal correlations between 
E and matter geometry. It is therefore interesting to also 
study these correlations. In Fig. 6, we show the correlations 
(cos[n(\&E — ^n)])i n = 1,2,3,4, as functions of the im- 
pact parameter at the four field points as before. There are a 
number of interesting features that differ from the magnetic 
field case: 1) there is a clear negative correlation (i.e. back-to- 
back) between "Je and v&i and it is most strong in the more 
central collisions and at the field point near the origin — this 
is understandable as the pole of with more matter will con- 
currently have more positive charges from protons generating 
E pointing in the opposite; 2) at field points away from the 
origin the ^e is strongly correlated to ^ 2 ( an d a ls° weakly 
correlated to $4) — the direction of E is more in-plane for the 
field point on x-axis while more out-of-plane for points away 
from x-axis; 3) there is also a weak correlation between ^e 
and ^3. Finally we have also checked the (e£') 2 -weighted 
correlations ((eE) 2 cos[n(* E - * n )])/((eE) 2 ) and find no 
visible correlation between the E magnitude and orientation. 



VI. DISCUSSION ON CHARGE-DEPENDENT 
MEASUREMENTS 

The studies of correlations between B and E fields and 
matter geometry suggest in general that there may be non- 
trivial charge distributions, particularly in azimuthal angles 
(e.g. charged dipole and quadrupole), induced by varied B- 
and E— induced effects. There are also other effects not re- 
lated the initial E and B fields, e.g., the local charge conser- 
vation effect [15] that can lead to nontrivial charge azimuthal 
correlations when coupled with various harmonic flows. It is 
therefore tempting to think about possible measurements that 
may fully extract information for the azimuthal charge distri- 
butions. In parallel to the measurements of (charge-inclusive) 
global azimuthal particle distributions that can be subse- 
quently Fourier-decomposed into various harmonic compo- 
nents, here we suggest a class of observables, the charged 
multiple vectors Q c n . Consider the measured charged hadrons 
in an event we can construct the charged multiple vector Q c n 



with magnitude Q c n and azimuthal angle ^ c n : 

QF n e* n * c »=Y i q i e* n * (6.1) 

i 

where the summation runs over all particles with qi and <fri 
the electric charge and azimuthal angle of the i-th particle. 
This idea generalizes the earlier proposal of the charged dipole 
vector analysis suggested in [19] in the context of CME ob- 
servables. One may introduce properly p t -weighed defini- 
tion. One may also think of sub-event version of this anal- 
ysis or possible multi-particle correlation improved version. 
We emphasize that these are different and independent from 
the existing Q n vectors related to flow measurements, i.e. 
Qn e m *™ = J^i e m< ^. The Q n is charge blind and includes 
all charges similarly while the Q c n takes the difference be- 
tween positive and negative charges therefore yields infor- 
mation on the charge distribution. The technical difficulty of 
Q c n measurements should be at the same level as previous Q n 
analysis and quite feasible, while clearly the Q c n analysis pro- 
vides orthogonal and unique information on the charge dis- 
tributions. With a joint Q n and Q c n analysis one can study 
the strength and azimuthal correlations among all harmonics 
and charged multipoles. Therefore it is very valuable to do 
a systematic charged multiple vector analysis, leading toward 
a quantitative "charge landscape survey" in heavy-ion colli- 
sions. 



VII. SUMMARY 

In summary, we have performed a detailed study of the 
event-by-event fluctuations of both the azimuthal orientation 
of the magnetic and electric fields as well as the matter geom- 
etry (which is specified by the participant planes of a series 
of harmonics) in the initial condition of heavy-ion collisions. 
Such fluctuations suppress the azimuthal correlations between 
the magnetic field and the second harmonic participant plane 
particularly in the very central and very peripheral collisions, 
while leaving a window around b = 8 ~ 10 fm with still 
sizable correlation between the two. We have further studied 
similar azimuthal correlations between B and other harmonic 
participant planes, and found similar yet weaker correlation 
with $4 while no correlation with ^1 and ^3. The correlation 
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between and ^4 may indicate that the CME can also con- 
tribute to the charge-pair correlation (cos[2((/) 1 + 4> 2 — 2^4)]) 
as opposed to the assumption in Ref. [35]. Examination of 
these correlations at different field points shows notable de- 
pendence on spatial positions. For completeness we have pre- 
sented similar studies for electric field E azimuthal correla- 
tions with matter geometry which show quite different pat- 
terns from that for the B field. We have evaluated the im- 
pact of such azimuthal fluctuations and correlations on a num- 
ber of observables related to magnetic-field induced effects 
in heavy-ion collisions: the charged pair azimuthal correla- 
tions, the charge dependent elliptic flow of pions, as well as 
the azimuthal anisotropy of soft photons due to conformal 
anomaly. Specifically we have quantified the modification 
factors Ri = (cos(24'b)) and the (eB) 2 -weighted reduction 
factor R 2 = ((eB) 2 cos(2* B )>/((eB) 2 ) at different central- 
izes (see Figs. 4 and 5) and found sizable reduction in both 
very central and very peripheral collisions. From these results 
we conclude that the optimal centrality class for observing the 
above mentioned B-induced effects corresponds to impact pa- 
rameter range b ~ 8 — 10 fm. The qualitative conclusion 
should hold also for Pb + Pb collisions at LHC. 

We end with a few pertinent remarks: 
(1) In our computations (as well as all the previous compu- 
tations [2-8]), the classical Lienard-Wiechert potentials have 
been used. One may worry about quantum corrections to 
the field equations as the magnitude of the electromagnetic 
field is much larger than the electron and light quark masses. 
In principle a calculation including all relevant QED pro- 



cesses is needed, but a one-loop Euler-Heisenberg effective 
lagrangian (see Ref. [36] for review) may give a good indica- 
tion on the size of such corrections. At strong-field limit the 
field equation derived from this lagrangian can be regarded 
as linear Maxwell equations with the renormalized charge 
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1 



In- 



An estimate implies 



that the computed electromagnetic field may be amended only 
by a few percent even for e|F| ~ IOOto 2 . 

(2) As shown in Refs. [5, 8], the magnetic field magnitude 
bears strong time dependence. Furthermore, the electromag- 
netic response of QGP could dramatically modify the time 
evolution of the B and E fields. It would be important to 
incorporate these into realistic modeling of effects induced by 
these fields, as recently attempted in [37] to incorporate the 
time evolution into the event-by-event estimates of the CME 
by using the parton-hadron-string-dynamics approach. 

(3) As suggested in Ref. [38], in order to disentangle the ef- 
fects of the magnetic field and the collective flow, it is very 
useful to study the U + U collisions where the geometry be- 
comes a nontrivial level arm. It will be interesting to study in 
the future the various correlations among the field strengths, 
field azimuthal orientations with the matter geometric eccen- 
tricities as well as the harmonic participant planes. 
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